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Abstract 

The magnitude of the CKM matrix element |Vub| is determined by measuring the in- 
clusive charmless semileptonic branching fraction of beauty hadrons at OPAL based on 
b — >■ Xu£z/ event topology and kinematics. This analysis uses OPAL data collected between 
1991 and 1995, which correspond to about four million hadronic Z decays. We measure 
Br(6 — * Xu^z/) to be (1.63 ± 0.53 llo]^^) ^ 10"'^. The first uncertainty is the statistical error 
and the second is the systematic error. From this analysis, |Vub| is determined to be: 



|Vub| = (4.00 ±0.65 (stat) 



:°:^^ (sys) ± 0.19 (HQE)) x 10^1 



The last error represents the theoretical uncertainties related to the extraction of |Vub| 
Br (6 — i> Xu^z/) using the Heavy Quark Expansion. 
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1 Introduction 



The CKM matrix [l|] describes the relation between quark weak and mass eigenstates, with 
the element Vub describing decays of the b to u quark. Its magnitude, |Vub|, can be cal- 
culated by measuring the inclusive b — *> u semileptonic decay rate. Given that the branch- 
ing fraction of inclusive b ^ u semileptonic decay is of order 10~^, a large number of b 
hadrons are required to measure |Vub|- The dominant background to 6 — >■ Xu^z/ comes from 
b -^ Xciu decays because the branching ratio of 6 — > X^^u is more than 50 times greater than 
that oi b —^ Xu^z/. Here the lepton i refers to either an electron or a muon, and b denotes all 
weakly decaying b hadrons^Q^. Xu and Xc represent hadronic states resulting from a b quark 
semileptonic decay to a u or c quark respectively. The determination of |Vub| depends on 
the b to u and b to c semileptonic decay models. 

The inclusive method developed by ARGUS and CLEO ^ is to extract |Vub|/|Vcb| 
from the excess of events in the 2.3 to 2.6 GeV/c region of the lepton momentum spectrum 
in the B meson rest frame, where the b -^ Xc^z/ contributions vanish. This technique uses 
only a small fraction of the lepton phase space and so has considerable model dependence 
in extrapolating to the entire lepton spectrum in the B rest frame. In addition, since the 
LEP experiments can not precisely determine the B meson rest frame, this method is not 
appropriate for the LEP experiments. Instead, at LEP, |Vub| or |Vub|/|Vcb| is extracted using 
a larger portion of the lepton spectrum as well as other kinematic variables. The inclusive 
measurement of the branching fraction of the b -^ Xu^z/ decay has been performed at LEP 
by ALEPH §, DELPHI § and L3 §. 

The theoretical uncertainty for the value of |Vub| extracted from a measurement of inclu- 
sive b — > Xu£z/ branching fraction differs from that extracted from measurements of exclusive 
b ^ u semileptonic decay rates. A recent theoretical study concludes that there is a 5% 
theoretical uncertainty on |Vub| values derived from b -^ X^lv inclusive measurements 0, 
using the Heavy Quark Expansion. There is a 15% theoretical uncertainty associated with 
I Vub I values extracted from measurements of the exclusive branching fractions B -^ Triu 
or B — i> piu 1^, interpreted within the framework of the Heavy Quark Effective Theory 
(HQET). 

In this paper, we describe the determination of |Vub| using the inclusive b — > X^^z/ decay 
rate from the OPAL data taken at center of mass energies near the Z resonance. The 
event preselection, the b — > X^iu decay models and the neural network used to separate 
b —>■ X^iu from the background will be discussed in detail in the following sections. 



2 The OPAL detector, data and Monte Carlo samples 

The OPAL detector is a multi-purpose Air spectrometer incorporating excellent charged and 
neutral particle detection capabilities. The OPAL detector is described in detail elsewhere . 
A brief description is given here. The central tracking system consists of a silicon microvertex 
detector, a vertex chamber, a jet chamber and z chambers. The momentum of tracks and 



^Charged conjugate states are implied if not stated otherwise. 



the primary and secondary vertex position are reconstructed by the central tracking system, 
which is located inside a solenoid. The solenoid provides a magnetic field of 0.435T. Outside 
the solenoid is the electromagnetic calorimeter, which is composed of lead glass blocks and is 
used to measure the energies and positions of electrons and photons. The hadron calorimeter 
lies outside the electromagnetic calorimeter and is used to measure the energy of hadrons 
emerging from the electromagnetic calorimeter and assists in the identification of muons. 
The outermost OPAL detector is the muon detector which consists of a system of barrel and 
endcap muon chambers. A large fraction of muons with momenta less than 2 GeV/c are 
absorbed by the other detectors or the iron shielding before reaching the muon chambers. 

The current analysis uses OPAL 1991 to 1995 data, collected near the Z resonance, com- 
prising about four million hadronic Z decays. Monte Carlo simulated events were generated 



using the JETSET 7.4 [jTO[ generator, with parameters described in [^. Approximately 
five million hadronic Z — ^ bb decays were generated to study the b -^ Xc^z/ decay and the 
h ^ c ^ i cascade decay. Six million hadronic Z — ;► qq (where q can be u, d, s, c and b) 
decays were generated to study the leptons from primary charm quarks and light quarks. 
Two hundred thousand events from a 6 — »• Xu^i^ hybrid model |]T2| were produced to simulate 



the b ^ u semileptonic decay. The hybrid model will be described in detail in Section 3.1. 



3 Signal and background simulation 

The b to u semileptonic decay and background simulation are described below. The b to u 
semileptonic decay and background simulated events are passed through the full OPAL de- 
tector simulation |13| to produce the corresponding response. For this paper, the production 
fractions of B"*", B'^, B^ and Ab in Z decays were adjusted to reproduce those given by the 
Particle Data Group |jl4 . 



3.1 The b -^ X^£i^ hybrid model 

Several theoretical models have been proposed for the b -^ y^Jlv decay. Exclusive bound- 



state models |T5|-[18| approximate the inclusive b -^ y^^tv lepton spectra by summing con- 
tributions from all the exclusive final states. These exclusive models do not include all the 
possible final states nor any non-resonant states and therefore yield an incomplete predic- 
tion of the inclusive lepton momentum distribution, especially in the region of high hadronic 
invariant mass. The inclusive free quark models |lT9|-p^ treat the heavy quark as a free 
quark and the final state as a quark plus gluons. Free quark models are known to give 
poor agreement with experiments at low u quark recoil momentum. Therefore, a hybrid 



model [|12| has been proposed to model the b -^ y^^lv decay by using the exclusive model in 



the lower hadronic invariant mass region and using the inclusive model in the higher hadronic 



invariant mass region. The ISGW2 model [0 is used as the exclusive part of the hybrid 



model. The ACCMM model 0, combined with the W decay model plus JETSET 
fragmentation, is used as the inclusive part of the hybrid model. Since the ISGW2 exclusive 
model includes the exclusive resonant final states IS, 2S and IP up to 1.5 GeV/c^ in the 
hadronic mass, the boundary between the inclusive and exclusive parts of the hybrid model 



is placed at the hadronic invariant mass of 1.5 GeV/c^. The relative normalization of the 
inclusive and exclusive parts of the hybrid model is determined by the inclusive model. This 
hybrid model is only applied to decays of B mesons. There are no theoretical predictions for 
b to u semileptonic transitions of b baryons. The exclusive transitions of the b baryons in 
the OPAL tune of JETSET [0,|n[ are used. 



In order to estimate systematic uncertainties due to modeling of the inclusive spectrum, 
alternative models are also studied. Signal events were generated with the QCD universal 
function [pO|-p2[1 and parton ||2^ models. The invariant mass distributions of the hadronic 
recoil uq system are shown in Figure |l]a for the QCD universal function, ACCMM and 
parton models. The invariant mass distribution of the hadronic recoil uq system for the 
hybrid model is shown in Figure |l]b. 
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Figure 1: a,b. The uq invariant mass distributions, a using the QCD universal function, 
ACCMM and parton inclusive models, b using the hybrid model. Only the portion of the 
uq invariant mass above 1.5 GeV/c^ from the inclusive model in a is used in the hybrid 
model. The boundary between the exclusive model (left arrow) and the inclusive model 
(right arrow) in the hybrid model is indicated by the dashed line in b. 



3.2 Background simulation 



The ACCMM model |T9| is used to describe the lepton spectrum of 6 — ;> Xciu and b — > 
c ^ i decays. The fragmentation function of Peterson et al. |^5|] is used to describe the b 
quark and c quark fragmentation. The branching fractions of B° — >■ T)~t^u, B° — > D*~£"'"i/, 
B"^ — > D'^f+z/, B+ — »• D*''£"'"z/ and Ab — > AcX^i/ were modified to reproduce those given by 
the Particle Data Group [0. The Ab lepton momentum spectrum corresponding to -56% 
polarization |^6| was used. 



4 Event preselection 



A hadronic event selection |^^ and detector performance requirements are applied to the 
data. The thrust polar angle^Q* |cos^| is required to be less than 0.9 to ensure that the 
events are well contained within the acceptance of the detector. The selected events must 
pass the b identification, the lepton selection and the b semileptonic decay selection. All 
these selections are described in detail in the following sections. After all these preselections, 
the b -^ Xu^z/ decay purity is 1.3% and the main background is from b -^ Xciu decays. 



4.1 b identification 



A neural network algorithm [^ based on charged particle vertex information is used to 
separate the b flavour events from the other flavour events in each hemisphere. If either 
hemisphere passes this neural network selection, the event is selected. After this neural 
network selection, the b purity is more than 91% and the b identification selection efficiency 
is approximately 30% per hemisphere from the Monte Carlo simulation in which a branching 
fraction of 1.0 x 10~^ for the b —>■ X^ii^ transition is incorporated. Both hemispheres are 
searched for electron and muon candidates after the b identification. 



4.2 Lepton selection 



Electrons are identified by a neural network ^E using the track and calorimeter informa- 
tion. The electron momentum is required to be greater than 2 GeV/c. Electrons from 
photon conversions, 7 — >■ e"'"e~, contribute a significant background to the prompt electron 



samples. Another neural network is used to reject this background |^. The photon con- 
version background is reduced by 94% after the photon conversion neural network selection, 
whilst retaining 98% of the selected prompt electrons. After all these requirements, the re- 
sulting electron efficiency is approximately 74% with a purity of 94% within the geometrical 
acceptance. 



Muons are identified using reconstructed track segments in the muon chambers [ p8[] . 
The muon momentum is required to be greater than 3 GeV/c. The reconstructed tracks 
in the central detector are extrapolated to the muon chambers to see if they match the 
track segments reconstructed in the external muon chambers. The measured energy loss 
dE/dx is also required to be consistent with the expected value for a muon. After all 
these requirements, the muon selection efficiency is approximately 90% and the muon purity 
approximately 93% within the geometrical acceptance. 

Electron and muon momenta transverse to the direction of the jet containing the lepton 
are required to be greater than 0.5 GeV/c in order to reject leptons from light quark decay. 
The lepton is included in the calculation of the jet direction. The jet finding is based on the 



cone algorithm |29 



^ A right handed coordinate system is used, with positive z along the e~ beam direction and x pointing 
toward the center of the LEP ring. The polar and azimuthal angles are denoted by 6 and (j), and the origin 
is taken to be the center of the detector. 



4.3 b semileptonic decay selection 



A neural network |]30| based on lepton information is used to separate the b hadron semilep- 
tonic decays, b — > X^^z/ and b —>■ X^^iu, from non b semileptonic decays. The distributions 
of the neural network output variable are shown in Figure |^. After this neural network b 
semileptonic decay selection, the b hadron semileptonic decay purity is 97% and the efficiency 
is 65% for this neural network; the c ^ i events, where c is a primary quark, and h ^ c ^ i 
events are suppressed. 
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Figure 2: a,b. The b hadron semileptonic decay neural network output distributions, a for 
the b semileptonic decays and the scaled background from the Monte Carlo simulated events, 
here the background indicates all events excluding b — >■ X^^z/ and b — >■ X^^z/ , b comparison 
between the OPAL data and the Monte Carlo simulated events. The selected region is shown 
by the arrow in b. 



5 b — ^ Xu^^' neural network 



Because of the dominant b -^ Xc£z/ background, it is difficult to enrich the sample in 
b -^ X^j^iu decays using only one kinematic variable. A multi-layered feed- forward artifi- 



cial neural network based on the JETNET 3.0 program [^ is used to enrich the sample in 
b — > Xu^z/ decays. There are four layers in this neural network. The neural network structure 
is 7-10-10-1. In the first layer, seven variables are used as inputs to the neural network. The 
last layer is the neural network output variable. A figure of merit [^] is used to determine 



the discrimination power of these seven variables in separating two classes of events, i.e. 
signal and background. The higher the figure of merit, the better the separation between 
the two classes. Over twenty kinematic variables were initially considered as inputs to the 
b —>■ Xu£z/ neural network. Only seven variables are selected as inputs to the b —>■ X^iu neural 
network based on good separation between b -^ Xu£z/ and background and good agreement 
between data and Monte Carlo simulated events. These seven input variables, in order of 
decreasing figure of merit, are: 



1. the invariant mass of the most energetic final state particle combined with the lepton, 

2. the lepton energy in the b hadron rest frame, where the b hadron energy and momentum 
are estimated using the techniques described in p3| , 



3. the lepton momentum transverse to the jet axis (the jet axis calculation includes the 
lepton), 

4. the transverse momentum of the most energetic final state hadron with respect to the 
lepton direction (assuming all hadrons are pions), 

5. the rapidity of the most energetic final state hadron calculated with respect to the 
lepton direction (assuming all hadrons are pions), 

6. the fraction of the reconstructed b hadron energy carried by the lepton, 

7. the reconstructed hadronic invariant mass, Mx, which is calculated by: 

Mx^ = E(WiEO' - E(WiPi)'' (1) 

i i 

where i denotes all hadronic tracks and clusters. Wi is the probability that the i**^ 
hadronic track or unassociated cluster comes from b decay and is calculated using 



the techniques described in ^^. Ei and pi are the energy and momentum of the i 



;th 



hadronic track or neutral cluster. 



Only the tracks and clusters from the same jet as the lepton are included in the calculation of 
these seven input variables. The seven input variable distributions for the b — »• X^^z/ and the 
background in the Monte Carlo simulation are shown in Figure |^. The agreement between 
the Monte Carlo simulated events and the OPAL data for these seven variables is shown in 
Figure ^. 

Twelve thousand b -^ X^^lu events, which were simulated with the hybrid model and 
have passed the event preselection, and the same number of background events from the 
multi-hadron Z — i> qq Monte Carlo simulation after the preselection are used to train the 
b -^ Xu£z/ neural network. Two other samples of signal and background events of the 
same size are used to test the neural network performance. The neural network output 
distributions from b — > X^^z/ and background are shown in Figure |^a. 

The background composition from the b — > X^^i/ neural network is shown in Figure |^. 
Ninety percent of the background in this analysis comes from the b -^ X^^z/ decay, 6.8% from 
the b — s> c decay with the c subsequently decaying to a lepton. Another 0.6% comes from 
the c —>■ i decay in which the c quark is the primary quark. Other background processes 
make up the remaining 2.6%, of which 36% is from the h —* r decay with the r subsequently 
decaying to an electron or a muon, and most of the rest of the background is from a pion or 
a kaon misidentified as an electron or a muon. 
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Figure 3: 1-7. Comparison between the signal b -^ X^iu and the background in the Monte 
Carlo simulation for the seven b — > X^^iu neural network input variables. The b —>■ Xu£z/ signal 
and background are normalized to unity. The input variables in plots 1 to 7 are in the same 
order as the input variables defined in the text of Section B. 
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Figure 4: 1-7. Comparison between the OPAL data and the Monte Carlo simulated events 
for the seven neural network input variables. A branching fraction of 1.63 x 10~^ for the 
h — >■ Xu£z/ transition is incorporated in the Monte Carlo simulation for comparison with the 
OPAL data. The input variables in plots 1 to 7 are in the same order as the input variables 
defined in the text of Section ^ 
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Figure 5: a,b. The b -^ y^^lv neural network output distributions, a for h 
background, b for different background components. 
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6 Extraction of Br(6 -^ X^iiy) 



The branching fraction of 6 — *> X^iu decay can be obtained from the best fit of the Monte 
Carlo simulated events to OPAL data based on the b -^ X^iu neural network output distri- 
butions. Br(6 -^ Xu^z/) is extracted from the b -^ X^iu neural network output distributions 
by minimizing: 



X 



E 



[N 



data 
k 



N. 



data 



'x f, 



MCk 



+ (1 - x)l 



MCb 



'■)? 



Tvrdata 



(2) 



where N|^^*^ is the number of events from the data in the k*'^ bin of the neural network 
output. Ndata is the total number of events in the data after preselection. The free parameter 
X is the fraction of signal events in the data after preselection, which can be converted to 
Br(6 -^ X^iu) based on the number of signal events and the number of background events 
in the Monte Carlo simulation after preselection, fj^ ''^ is the fraction of simulated signal 
events in the k**^ bin of the b -^ X^lv neural network output with respect to the total number 
of simulated signal events after preselection, fj^ ^'^ is the fraction of simulated background 
events in the k**^ bin of the b -^ XJlv neural network output with respect to the total number 
of simulated background events after preselection. Here the background includes b -^ Xc£z/, 
b^c^£, c^£ and other contributions. The sum over the index k is performed from 
the neural network cut to the last bin in the neural network output distribution. The 
Br (6 — > X^iv) from the fit result x, as well as its statistical and systematic errors, depends 
on the b -^ XJlv neural network cut. The resulting Br(6 -^ X^lv) is stable, with variations 
less than 0.1 x 10~^, as the neural network cut varies in value from 0.3 to 0.7. A neural 
network cut of 0.7 is chosen to minimize the total relative errors and yields 



Br (6 -^ X, 



;i.63±0.53) X 10" 



where the uncertainty is the statistical error only. 



In Figure |^a, the neural network output distribution from data and the Monte Carlo 
simulation events with no b ^ Xu^z/ semileptonic decay is shown and the excess of events in 
the data can be seen in the last bin. Here the distribution of Monte Carlo simulated events 
with no 6 — >■ Xu^z/ transitions is normalized to the same number of entries as the data. The 
X^/ndf is 14.6/9, which corresponds to a 10% confidence level, assuming no contributions 
from b -^ X^iu transition. Here the x^ is calculated by summing over all bins in the neural 
network output distribution. When the b — > Xu£z/ transition is incorporated in the Monte 
Carlo simulation with a branching fraction of 1.63 x 10~^, the Monte Carlo simulation 
agrees much better with the data, as can be seen in Figure ^. The x^/ndf is then 8.3/8, 
corresponding to a 41% confidence level. 

The data after subtracting the background from the Monte Carlo simulated events agree 
well with the simulated b -^ X^iu signal within statistical errors, which is shown in Figure ^ 
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Figure 6: a,b. The neural network output distributions for data and Monte Carlo simulated 
events, a with no 6 — *> Xu£z/ transition in the Monte Carlo simulated events, b with a 
branching fraction of 1.63 x 10~^ b — >■ Xu^i/ decay incorporated. The distribution of Monte 
Carlo simulated events is normalized to the data for both plots. 

A series of cross checks, dividing the lepton samples into electron and muon samples and 
dividing the data into two samples for the years 1991 to 1993 and the years 1994 to 1995, 
are performed. The Br(6 —>■ X^iu) results are consistent within statistical errors for these 
cross checks. 



7 Systematic errors 



The list of systematic errors is given in Table |1|. Unless otherwise specified, the systematic 
errors are estimated by varying each parameter described by one standard deviation and 
taking the corresponding largest errors. The resulting systematic errors in Table ^ are 
discussed in detail: 



Error Source 


Variation or 


ABr(6 ^ XJu) 




value and variation 


10-3 


Fragmentation (xE)b 


0.702 ± 0.008 @] 


+0.28 
-0.32 


Lepton spectrum (b -^ c) 


ISGW** 36 , ISGW 16 


+0.18 
-0.29 


MC statistics 


(see text) 


±0.22 


b and c hadron semileptonic decay 


(see text) 


±0.19 


MC modeling 


(see text) 


±0.19 


b —>■ Xu£z/ modeling error (hybrid) 


(see text) 


±0.19 


b -^ X^iu modeling error (inclusive) 


Parton 23 , QCD |0| 


±0.14 


b -^ X^iu modeling error (exclusive) 


ISGW2 18 , JETSET 


±0.07 


Tracking resolution 


±10% 2S| 


±0.07 


c hadron decay multiplicity 


(see text) 


±0.07 


Ab production rate 


(11.6 ± 2.0)% ll 


^0.04 


Ab polarization 


(-561^?)% 1 


±0.03 


Electron ID efficiency 


±4% H 


^0.04 


Muon ID efficiency 


±2% ig 


^0.03 


Electron fake rate 


±21% 2| 


^0.02 


Muon fake rate 


±8% H 


^0.01 


Br(6 ^ Xrz/,) 


(2.6 ± 0.4)% 14 


±0.01 


b lifetime 


(1.564 ± 0.014) ps U 


< 0.01 


Rb 


0.2178 ± 0.0017 


< 0.01 


Total 




+0.55 
-0.62 



Table 1: Systematic errors for Br (6 —>■ Xu 
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Figure 7: a,b. The neural network output distributions, a Data after subtracting the back- 
ground from the Monte Carlo simulated events (points) show agreement with the simulated 
b -^ Xu-^z/ signal (solid histogram), b Data after subtracting the Monte Carlo simulated 
events with a branching fraction of 1.63 x 10^'^ b —>■ Xy.iu decay incorporated. Here the 
error bars include the statistical error from data and Monte Carlo simulated events. 



b quark fragmentation: Many parameterizations have been suggested to describe the 
heavy quark fragmentation process. The Peterson function EH] is used here to simulate 
the b and c fragmentation in the Monte Carlo simulation. The systematic error in the 
b quark fragmentation is estimated by varying the b hadron mean scaled energy (xE)b 
within the experimental range 0.702 ± 0.008 recommended by the LEP Electroweak 
Working Group [^. This value is consistent with a recent determination of (a;E)b = 



0.714 ± 0.009 from SLD [^. The systematic error is also estimated from the Collins 



and Spiller fragmentation function |38] and Kartvelishvili fragmentation function [39 



The uncertainties in c quark fragmentation can be neglected because the background 
from c —>■ i, where c is a primary quark, is very small. 

Xc^i^ lepton momentum spectrum modeling: Different decay models are used to 
predict the lepton spectrum in the b hadron rest frame for the b — > X^iu decay. Al- 
though all models are derived for B'' and B+ semileptonic decay only, they are extrapo- 
lated to the Eg and Ab semileptonic decay. This will be correct in the simple spectator 
model and is a reasonable approximation for this analysis. The lepton spectrum from 
the ACCMM model [^ is used as a base model for the b — > Xc^u decay. The sys- 
tematic errors due to 6 ^ X^iiy lepton momentum spectrum modeling are estimated 
from the ISGW ^^ and ISGW** EE models as prescribed by the LEP Electroweak 



Working Group |35 



The lepton spectrum from the b — »• c — ;► £ decay in the ACCMM model is different 
from the lepton spectrum in the ISGW model. The systematic error due to the shape 
of the b ^ c — > £ lepton spectrum is calculated and is found to be negligible. The 
lepton spectrum from the c — *> £ decay is varied from the ACCMM model to the ISGW 
model, where the c quark is a primary quark from Z decay. The systematic error is 
calculated and found to be neghgible. 



Monte Carlo statistics: The systematic uncertainty due to tlie limited Monte Carlo statis- 
tics is ±0.22x10-3. 

b and c hadron semileptonic decay branching fractions: The systematic error is es- 
timated from the uncertainties of the branching fractions of B ^ D£u, B -^ D*£z/, 
B — > D**£i/ and Ab -^ AcX£z/. There is a 6.8% background contribution from the 
h -^ c ^ i decays and a 0.6% background contribution from the c — >■ £ decays. The 
systematic error is also estimated from the uncertainties of the branching fractions of 
the h -^ c ^ i decays. A summary of these systematic errors from the uncertainties 
of b hadron and c hadron semileptonic decay branching ratios is shown in Table ^ 
The Br(B — > D**^z/) in Table || is obtained by averaging the Br(B -^ D**£i/) from 
ARGUS gg], ALEPH gT|], DELPHI [g^ and the total B semileptonic decay branching 
fraction subtracting the contribution from B to D and D* semileptonic decay, de- 
scribed by the LEP, CDF and SLD Heavy Flavour Working Group ^^. For the decay 
of B ^ D**iu, in which D** refers to Di, D2, D2 and Dl, the branching ratio for each 
specific D** final state is not well measured. For this analysis, the narrow final states 
of D** in B — > D**£z/ are replaced by the broad states and then vice- versa to check the 
sensitivity of the Br (6 -^ Xu£i/) to the relative ratio of the narrow and broad states of 
D** in B — > D**iu. The effect on the Br(6 -^ Xu£z/) is found to be negligible. 



Error Source 


Variation 


ABr(6^Xu£z/)(10-3) 


Br(B° -^ D-£+u) 


(2.10 ± 0.19)% g 


^0.02 


Br(BO -^ D*-i+u) 


(4.60 ± 0.27)% ig 


±0.03 


Br(B+ ^ T)H+u) 


(2.15 ± 0.22)% g 


^0.06 


Br(B+ ^ D*0£+z/) 


(5.3 ± 0.8)% 14 


±0.04 


Br(B -^ D**£u) 


(3.04 ± 0.44)% g 


±0.16 


Br(b ^c^i) 


(8.4 tiii)% m 


^0.02 


Br(Ab -^ A^Xiu) 


(7.9 ±1.9)% in 


TO.06 


Total 




±0.19 



Table 2: Systematic errors for Br(6 — 
semileptonic decay branching ratios. 



X„ 



from uncertainties of the b hadron and c hadron 



Monte Carlo modeling errors: The systematic error for the Monte Carlo modeling er- 
rors is estimated by re-weighting each input variable distribution in the Monte Carlo 
simulation to agree with the corresponding data distributions. A branching fraction 
of 1.63 X 10~3 for the b -^ X^iu transition is incorporated in the Monte Carlo sim- 
ulation as shown in Figure ^. This gives a conservative estimation of the systematic 
uncertainty due to the modeling of the input variables. 

b — > Xu^i^ modeling error from the hybrid model: The boundary between the inclu- 
sive and exclusive regions in the hybrid model is varied from 1.5 GeV/c^ to 0.9 GeV/c^. 
This conservatively estimates the systematic error arising from the placement of the 
boundary between the inclusive and exclusive models. This produces a uncertainty of 
±0.19x10-3 for Br(6 -^ XJp). 

b — > Xu-^i^ inclusive model: The ACCMM model is the base inclusive model. The QCD 
universal function model and the parton model are used to estimate the systematic 



errors in the inclusive part of the b -^ X^^iu hybrid model. This gives a change of 
-0.14 X 10~^ for the QCD universal function model and +0.02 x 10~^ for the parton 
model for the branching ratio of 6 — >■ X^iu. The largest variation of Br(6 — > X^iu) 
from these models is taken as the systematic uncertainty of Br (6 -^ Xu^z/) from the 
inclusive model. 

Xu-^i^ exclusive model: The 1SGW2 model is the base exclusive model. The model 
implemented in the OPAL tune of JETSET |TI| Monte Carlo simulation, which has 



the u quark and the spectator quark forming one single hadron in the final state, is 
used to estimate the systematic error in the exclusive part of the b -^ X^iu hybrid 
model. 

Tracking resolution: The systematic error due to the uncertainties of the detector res- 
olution is estimated by applying a ±10% variation to the r-0 track parameters and 
an independent ±10% variation to the analogous parameters in the r-z plane to the 
Monte Carlo simulated events Ej 



c hadron decay multiplicity: The systematic error of the Br(6 -^ Xu^z/) associated with 
the c hadron decay charge multiplicity is estimated using the average charged track 
multiplicity of D"*", D°, D^ decays as measured by MARK III |43]. The systematic 



uncertainty of the Br(6 — > Xu£z/) is ±0.07x10 ^ from the uncertainty of c hadron 
decay multiplicity. 

Ab production rate: The PDG ||14| gives the production fraction of B"*", B*^, B^ and Ab in 
Z decay as (38.9 ± 1.3)%, (38.9 ± 1.3)%, (10.7 ± 1.4)% and (11.6 ± 2.0)% respectively. 
The neural network output variable distributions among B+, B° and B^ are similar and 
the systematic effects caused by the uncertainties of the production fraction of B+, B° 
and Bg are negligible. Due to the difference of the neural network output variable 
distributions between Ab and B mesons, the fraction of Ab is varied by one standard 
deviation to determine the corresponding systematic error. 

Ab polarization: The lepton momentum spectrum from Ab semileptonic decays depends 
on the degree of Ab polarization. The systematic uncertainties are estimated by using 
the Ab polarization ranging from -13% to -87%, as 95% [^ confidence level limits, 
which are converted to one standard deviation in Table |T]. 

Lepton identification efficiency: The number of selected events in the signal and back- 
ground depends on the electron identification efficiency and the muon identification 
efficiency. The electron identification efficiency has been studied using control samples 
of electrons from e+e~ — > e+e" events and photon conversions, and is modeled to a 
precision of 4% [Q. The muon identification efficiency has been studied by using 



the muon pairs produced in two-photon collisions and Z -^ //"'"/i events, giving an 
uncertainty of 2% |3^ . 



Lepton fake rate: Fake electrons in the electron sample are primarily from charged hadrons 
(mainly charged pions) misidentified as electrons and from untagged photon conver- 



sions. The uncertainty associated with electron misidentification is ±21% p8|. The 
muon fake rate is studied from K^ -^ tt+tt^ and t ^ Stt decay. The uncertainty of the 
fake muon rate is estimated to be ±8%. 



b — y XTi>T branching ratio: One important composition in the "other" background in 
Figure |p3 results from a b quark semileptonic decay to a r lepton, with the r lepton 
subsequently decaying to an electron or a muon. The branching ratio of 6 — >■ Xru^ 
is (2.6 ± 0.4)% |T^. The systematic error is estimated using the uncertainties of the 
h —>■ yLrUr branching ratio. 

Uncertainty of the b lifetime The average b hadron lifetime is measured to be (1.564 
± 0.014) ps |T^. The uncertainty in b lifetime results in a negligible uncertainty in 
Bi (b^XJu). 

Uncertainty of Rb: The fraction of Z — > bb events in hadronic Z decays, Rb, is measured 
to be 0.2178 ± 0.0017 [l3|. The uncertainty in Rb results in a negligible uncertainty 
in the background composition. 



8 Conclusion 

The branching fraction of the inclusive b -^ X^iu decay is measured to be: 

Br(6 -^ XJu) = (1.63 ± 0.53 (stat) +°;| (sys)) x 10^1 

The first error 0.53 is the statistical error from the data only. The errors associated with the 
limited statistics of the Monte Carlo sample are included in the systematic error. This result 
is consistent with similar measurements from ALEPH, DELPHI and L3, the other three LEP 
experiments, shown in Table |^. In Table |^, the first error in Br(6 —>■ Xu£i/) combines the 
statistical error from the data and limited Monte Carlo statistics as well as the uncorrelated 
systematic uncertainties due to experimental systematic errors, such as detector resolution 
and lepton identification efficiency. The second error contains the systematic uncertainties 
from the b -^ Xciu Monte Carlo simulation models. The third error contains the systematic 
uncertainties from the b — > X^iu models. The second and third errors are correlated between 
the various experiments. 



Experiment 


Br(6 -^ XJu) (10-3) 


Ref 


ALEPH 


1.73 ± 0.56 (stat+det) ± 0.51 (b ^ c) ± 0.2 (b -^ u) 


i 


DELPHI 


1.69 ± 0.53 (stat+det) ± 0.45 (b ^ c) ± 0.2 (b -^ u) 


i 


L3 


3.3 ± 1.3 (stat+det) ± 1.4 (b ^ c) ± 0.5 (b -^ u) 


i 


OPAL (This analysis) 


1.63 ± 0.57 (stat+det) tHi (b -^ c) ± 0.25 (b -^ u) 





Table 3: Br(6 -^ XJv) resuhs from ALEPH, DELPHI, L3 and this analysis. 



|Vub| can be obtained from Br(6 -^ X^dv) ^,^ with inputs slightly revised as described by 
the LEP Heavy Flavour Working Group |^3[ in the context of the Heavy Quark Expansion |0|: 



|Vub| = 0.00445 X 



'Br(&^Xu£z/)1.55psV 



0.002 



Th 



1 x(l±0.010pert±0.030i/„.3±0.035^J, (3) 



where the average b hadron hfetime Tb is equal to (1.564 ± 0.014) ps [O] from the LEP 
experiments. Thus, |Vub| obtained from this analysis is: 



|Vub| = (4.00 ±0.65 (stat) +°:^^ (sys) ±0.19 (HQE)) x 10"^ 

where the systematic error includes the b to u and b to c semileptonic decay modeling error, 
and the HQE error is purely the theoretical error from the Heavy Quark Expansion. This 
result is consistent with the |Vub| value from the CLEO exclusive measurement of (3.3 ± 0.8 
(total)) X 10-3 [|g. 
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